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^! Abstract 

We report measurements of the dependence on magnetic field and array 
size of the resistive transition of Josephson junction arrays with long-range 
interaction. Because every wire in these arrays has a large number of nearest- 

o : 

neighbors (9 or 18 in our case), a mean- field theory should provide an excellent 
description of this system. Our data agree well with this mean-field calcula- 
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tion, which predicts that T c (the temperature below which the array exhibits 
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■ macroscopic phase coherence) shows very strong commensurability effects and 
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scales with array size. 
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We report an experimental investigation of ordered Josephson junction Arrays with Long- 
Range Interaction (ALRI), of the sort originally proposed in the disordered limit by Vinokur 
et a/.0 Although such arrays had been fabricated by Sohn et a/J, the samples used in the 
present work for the first time have low enough critical currents and hence low enough 
screening to be in the regime well described by existing theoretical modelsi'S. 

These arrays consist of two perpendicular sets of N parallel superconducting wires, cou- 
pled by Josephson junctions at every point of crossing (see Fig. p. In this geometry, any 
horizontal (vertical) wire is nearest neighbor to all vertical (horizontal) wires, and next- 
nearest neighbor to all other horizontal (vertical) wires. Hence we term the interaction 
long-range. The number of nearest neighbors in these arrays is equal to the array size N. 
This is in sharp contrast to standard (short-range interaction) 2D arrays where the number 
of nearest neighbors (typically 4 or 6) is independent of array size. 

Arrays with long-range interaction were first proposed as a model for a spin glass for 
the case where the wires are positionally disordered and a sufficiently strong perpendicular 
magnetic field is appliedB. More recently, Chandra et aZ.i have shown that even for an 
ordered array, glassy behavior is expected in a very weak field (less than one flux quantum 
per row). The equivalent of "spins" in these ALRI are the phases of the superconducting 
wires. Above a transition temperature T c , the phases of the wires are uncorrelated. However 
when the array is cooled below T c , a transition to a macroscopically phase coherent state is 
predicted to occur. 

For an ordered array with long-range interaction in the limit of negligible screening, 
Sohn et alE have performed a mean-field analysis and computed the transition temperature 
T C MF (/) as a function of the applied field and array size. Because each wire has a large 
number of nearest-neighbors, a mean-field theory using the phase of each wire as a classical 
thermodynamic variable should provide a good description of this system. At zero field 
they find T c = NEj(T = T c )/2fce, where Ej(T) = hi c (T)/2e and i c (T) is the (unfluctuated) 
critical current of a single junction at temperature T. Note the unusual result that T c should 
scale with the size of the array. To keep T c of the array well below T™ ire , one requires 
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Nhi° c < Aek B Tf re , where i° c is i c (T = 0). The computations of Refs. [1], [3], and [4] only 
hold in the limit of negligible screening, where the array as a whole screens much less than 
one flux quantum <3> , and in the limit when phase gradients along the wires due to current 
flow are much less than phase drops at the junctions. The former condition can be written 
as N 2 L g i Q c <C $o, where L g is the geometric inductance of a cell in the array. The latter 
condition can be expressed as i^ nctlon <^ i™ ire . These three inequalities place very strong 
limits on the magnitude of i° c for given N . 

In the experimental work of Sohn and co-workers, Nhi° c / AchsT™ 1 ™ « 300 and 
N 2 L g i° c /<&o rs 10 3 . Hence their samples were not in the regime defined by the above men- 
tioned theories. We present here the first measurements of ALRI with critical currents small 
enough (of order 5 nA) to be in the limit of extremely weak screening, and to have an array 
T c well below the wire critical temperature. Our data show impressive agreement with the 
mean-field theory, including extremely strong commensurability effects. 

The samples consist of 0.25 /im wide Al wires {Tf re m 1.7 K) connected by Al-AlO^-Al 
junctions, fabricated as follows. A grid-like pattern of lattice constant 2 /xm is defined using 
electron-beam lithography on a PMMA-coated Si wafer. A three-angle shadow evaporation 
technique is used to deposit both sets of wires sequentially without breaking vacuum, using 
only the single lithography step. The evaporations are done at 45° to the substrate surface, 
but at different orientations with respect to the patterned channels. 30 nm of 99.999% pure 
Al are evaporated in the direction of one set of wires (the "horizontal" set). Al accumulates 
on the substrate only along those horizontal wires because the PMMA shadows the "vertical" 
wires. 150 mTorr of 2 is bled into the chamber, and an oxygen plasma is ignited for 
20 minutes to grow an AlO^ barrier. After pumping out the O2, the sample is rotated 
so that the second and third evaporations (30 nm of Al each) are done in the direction 
of the "vertical" wires, going "up" for the second evaporation and "down" for the third, 
to ensure that the vertical lines are continuous where they "climb" over the horizontal 
wires. A liftoff completes the process. This shadow evaporation technique yields very high 
quality underdamped junctions which are a major improvement over those from the previous 
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fabrication technique. □ 

The typical single junction resistance is Rff =70 kf2 which corresponds! to an unfluctu- 
ated critical current i° c of 5.6 nA, or Ej(T — 0) / ks — 0.13 K. Junction uniformity, measured 
from single junctions co-fabricated with the arrays, is approximately ±15%. The lead con- 
figuration is shown in the inset of Fig. |2|. Current is injected in the first wire of one set, and 
extracted from the last wire of that same set. We report data on two arrays: one consisting 
of 9 x 9 wires (8x8 cells) and the other of 18 x 18 wires (17 x 17 cells). 

The arrays are cooled to 315 mK in a 3 He cryostat within a double /i-metal shield 
which reduces the stray field to less than 50 mG. (A field of 5.2 G corresponds to / = 
$ ' ce ii /$ = 1 f° r our 2 jum spacing). Temperature stability is better than 3 mK below 2 K. 
A small magnetic field is applied perpendicular to the array using a solenoid surrounding 
the vacuum can of the cryostat. Screening by the array can be neglected because i° c is so 
small. Quantitatively, the ratio of the maximum flux screened by the array to the flux 
quantum is much less than one: N 2 L g i° c /<&Q «3 x 10 -3 -C 1, for N — 18 and where L g « 4 
pH is the geometric inductance of a single cell in the array, modeled as a superconducting 
square washer.! Considerable care was taken to ensure that the arrays are well shielded from 
RF and microwave radiation by the use of a shielded room, room-temperature low-pass LC 
Il-filters, cold resistors, and cold microwave filters.@ 

The current-voltage (I-V) curves for single junctions co-fabricated with the arrays do not 
show a well-defined critical current at 0.3 K because Ej < ksT, and hence a finite resistance 
is observed for all bias currents. The arrays on the other hand, consisting of many junctions 
in parallel, do show, at least at the lowest temperatures, a well-defined critical current and 
strong hysteresis, as expected from underdamped SIS junctions. Fig. ^| shows an I-V curve 
for the 17 x 17 array at / = 0. The two jumps in voltage to 2A/e and 4A/e (where A is the 
superconducting gap) correspond to all the junctions connected to one, then the other, of 
the current injection wires going normal. The unfluctuated zero-temperature array critical 
currents I® = Ni° c (~60 nA for the 9x9 wire array, ~100 nA for the 18 x 18 wire array) are 
so small that the measured J c will be significantly less than J° due to thermal fluctuations. 
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The measured I c actually corresponds to a jump from a finite-voltage (of order 1 /xV) phase- 
diffusion branch! to 2 A /e at a current value which is affected by damping as well as Ej and 
T. 

We therefore focus instead on the differential resistance Rd = dV/dl (at zero dc bias) as 
a function of field and temperature since it should be a better measure of the zero-current 
phase coupling of the array (and hence T c ). Rd is measured using a PAR 124 lock- in amplifier 
at 15.6 Hz with an excitation current of 0.2 nA (corresponding to ~ J c /10). Fig. ||| shows 
Rd vs. / = $ ce ^/$o plots for several temperatures from 0.3 K to 1.6 K for the 17 x 17 cell 
sample. The curves are not offset. Because Rd{f) is periodic in / and symmetric around 
/ = 1/2, we only plot Rd{f) for / ranging from to 1/2. Rd{f) displays minima at all 
commensurate fields where / = p/q, p and q being integers smaller than N. The / = 1/17 
and / = 1/16 states are not clearly resolved but all other commensurate states are clearly 
present (such as, for instance, all other multiples of 1/17, like 2/17 and 3/17). All the 
measured positions of the resistance minima are within less than 10~ 4 x $ from the ideal 
computed commensurate field values. We observe very similar behavior for the 8x8 cell 
array, with resistance minima at / = 1/8, 1/7, ... It is a characteristic feature of ALRI that 
such strong and detailed structure in the Rd{f) curve is visible. Standard 2D arrays and 
wire networks do not exhibit such richness of structure because they do not have the long- 
range order needed to support a stable vortex superlattice with such a large lattice constant 
(e.g. 17 cells). 

The deepest resistance minima occur at the most strongly commensurate states: / = 0, 
1/2, 1/3, 1/4. The shape of the Rd(f) curve is very similar near all of these states (see Fig. 
|3[). The full widths of the resistance dips (i.e. the field intervals between local maxima on 
either side of the dips) scale as 1/q, with q — 1, 2, 3 . . . Near integer / (e.g. / = 0, 1, 2 . . .) 
where q = 1, the resistance increases smoothly from f — n until / = n ± [l/(A r — 1)], i.e. 
the first adjacent commensurate state, giving a modulation-free half-width of 1/(N — 1). 
Corresponding behavior occurs near other strongly commensurate states. 

Rd(f) was measured for 20 temperatures between 0.315 K and 1.8 K, of which 11 are 
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shown in Fig. |3|. As the temperature is increased, Rd{f) increases and the relative am- 
plitudes of the resistance oscillations decrease until at higher temperatures (but with the 
wires still superconducting) Rd{f) saturates at the normal state resistance of the array 
R a ^ ray = 2R^f /18. In order to extract T c (/) from the Rd{f) vs. T data, we make use of 
the finite width of the resistive transition. We define the experimental T c using a resistive 
criterion; for each field value, Rd is plotted vs. T, and T c is taken to be the temperature at 
which Rd interpolates to eRjy, and e is a number between and 1. Automating this process 
produces the top two curves of Fig. [| of T c (f) for e = 0.5 (top curve) and e = 0.375 (middle 
curve). For values of e between approximately 0.4 and 0.8, the inferred T c scales almost 
linearly with e. 

The bottom curve is the result of a mean-field calculation of T , c Mi? (/) for a 17 x 17 cell 
array. T C MF (/) is the temperature above which the order parameter rji = (e 1 ^) is equal to 0, 
where <fii is the phase of the i th wire and the brackets denote a thermal average. There are 
no free parameters in this calculation, which consists of using an efficient scheme to find the 
largest eigenvalue of a 17 x 17 matrix given by Eq. (19) of Ref. [3] for each of one thousand 
field values shown. The eigenvalue problem is solved assuming a temperature-independent 
Ej, and T c Mi? is finally corrected to account for Ej(T), which varies by ~ 30 % over the 
temperature range of interest. 

The data and mean-field theory curves are in good agreement, both for the 17 x 17 cell 
array in Fig. [|, and for the 8x8 cell array (not shown). The maxima in the experimental 
T c (f) obviously occur at commensurate fields to the same high accuracy as the minima in 
the resistive data do, since the critical temperature was extracted from the Rd(f) curves. 
The mean-field theory also predicts local maxima in T C MF at all commensurate fields: we 
find that the positions of the clearly discernible maxima in the experimental T c (f) and T C MF 
agree to better than one part in 10 . As can be seen in Fig. ^, the lower resistance criterion 
gives better quantitative agreement with the mean-field theory T C MF , which is always below 
the experimental T c (and is defined slightly differently). We cannot use a resistive criterion 
of less than e ~ 0.37 over the whole field range because at low temperatures, for / ^ 0, 
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the array resistance saturates at a non-zero value (up to m 3 kf2 for incommensurate /), 
probably due to macroscopic quantum tunneling of the phases. 

In order to compare the experimental T c 's of the 18 x 18 wire and 9x9 wire arrays, we 
must first account for the temperature dependence of Ej in order to obtain T*, the transition 
temperature one would observe if Ej were constant and the same for both arrays. For / = 
we then obtain T*i 8 /T* 9 = 1.9, using e = 0.5 to determine T c for both arrays. This is very 
close to the theoretical value of 18/9 = 2, indicating that T c does indeed scale with system 
size N. 

In the absence of screening we can write the following simple expressions for the phases 
of each wire at x — ja and y = ka at zero applied current: 

V% = <p$ + 2irfNk 

<Pj = + 2irfkj. 

(Ph is the phase of the k th horizontal wire (constant along the wire) and ipj is the phase of the 
jth ver ^i ca i w j re (depends linearly on the position y along the wire), a is the lattice constant 
and the vector potential has been chosen as A = fx& /a 2 y . The only free parameter is 
A(Pq V = tfQ — <£q . The system energy E is 

JV-l 

E = - £ cos(v^-v^) 

k,j=Q 

The ground state energy is found by minimizing E numerically as a function of Aip^ v for 
each field. Once A<y?^ y is found, all the phase differences are then determined. The local 
extrema of both —E min (f) and T C MF (/) occur at exactly the same fields, with very similar 
relative amplitudes, indicating that the above simple expressions for the phases of the wires 
do indeed describe the phases very accurately. 

It is very difficult to probe the glassiness of this system using transport measurements 
because the phases unlock as soon as a small transport current is applied. Even though 
the arrays are biased well below I c , a finite voltage develops across the system because of 
phase diffusion. Phase diffusion is unavoidable in the small (i.e. low-capacitance) and weak 
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junctions required to conform to the model conditions. Since the phases are evolving as 
(dtp/dt) = 2eV dc /h, they cannot lock. Hence individual metastable states, the presence of 
which would confirm the presence of a glass, cannot be probed using our transport tech- 
nique. For instance, we observe the same I c at every field cool, while trapping into different 
metastable states should give a range of measured critical currents. Similarly, the T c we 
measure reflects an average over many configurations and thus reveals very little about the 
glassiness of the array. 

In conclusion, we have fabricated Josephson junction arrays with long-range interaction 
and extremely weak critical currents. A mean-field theory provides an excellent description 
of this system because every wire has a large number of nearest-neighbors (9 or 18 for the 
arrays presented here). Our data for T c (f,N) are in very good agreement with the mean- 
field calculation: we find that T c (f = 0) scales with system size and observe very strong 
commensurability effects. The array differential resistance at zero dc bias exhibits minima at 
all commensurate fields, displaying far more complex, but well understood, structure than 
standard 2D arrays or wire networks. 
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measurements and analysis, and J.M. Hergenrother, D. Davidovic, and M.S. Rzchowski for 
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JSEP Grant No. N00014-89-J-1023, and NSF Grant No. DMR-92-07956. 
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FIGURES 

FIG. 1. Schematic drawing of a 2 wire by 3 wire array with long-range interaction. There are 
Josephson junctions at every crossing point of the superconducting wires. 

FIG. 2. Voltage-Current plot at 0.315 K in zero field of the 17 x 17 cell array. The dashed 
line corresponds to sweeping current up, the solid line to sweeping current down. There is a finite 
slope at zero bias, too small to be seen on the graph. The inset is a schematic diagram of the lead 
configuration used for current injection and voltage measurement. 

FIG. 3. Plot of zero-bias differential resistance of the 17 x 17 cell array vs. normalized flux /, 
measured with a 0.2 nA ac excitation, for selected temperatures. The curves are not offset. From 
the lowest to the highest curve, the temperatures are: 0.417 K, 0.702 K, 0.797 K, 0.959 K, 1.047 
K, 1.13 K, 1.212 K, 1.288 K, 1.39 K, 1.51 K, and 1.69 K. The local minima in resistance occur 
within 10~ 4 x $o of all the commensurate flux values, i.e. at all f = p/q where p and q are integers 
between 1 and 17. 

FIG. 4. Plot of the temperature T c corresponding to the onset of macroscopic phase coherence 
vs. normalized flux / for the 17 x 17 cell array. The top two curves (data) are computed from the 
differential resistance vs. field data by using a resistive transition criterion for T c of 0.5i?^ ray (top 
curve) and 0.375i?^ raj/ (middle curve). The lower curve is the result of a mean-field calculation of 
T c MF (f). 
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